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In this paper, we study boundary-induced phase transitions in a particle non-conserving asym- 
metric simple exclusion process with open boundaries. Using boundary layer analysis, we show that 
the key signatures of various bulk phase transitions are present in the boundary layers of the density 
profiles. In addition, we also find possibilities of surface transitions in the low- and high- density 
phases. The surface transition in the low-density phase provides a more complete description of the 
non-equilibrium critical point found in this system. 



I. INTRODUCTION 

Studies on a certain class of non-equilibrium systems, 
namely driven diffusive systems, have revealed many fea- 
tures that are unexpected in systems in thermal equilib- 
rium P, For example, it is well known that although 
in thermal equilibrium, one dimensional systems with 
short-range interaction cannot exhibit phase transition or 
spontaneous symmetry breaking, certain driven diffusive 
systems do exhibit such phenomena even in one dimen- 
sion Asymmetric simple exclusion processes (ASEP), 
which involve biased hopping of particles in one direction 
with hard-core exclusion along a one-dimensional lattice, 
fall in the class of driven diffusive systems and despite 
its simplicity, such processes are capable of exhibiting 
variety of phenomena 0, 01 E3 Violation of detailed bal- 
ance due to the presence of a finite particle current makes 
such systems so special. This, for example, is reflected 
in the boundary induced phase transitions in ASEP on 
one-dimensional finite lattice 0, where, unlike equilib- 
rium, the effect of boundaries propagates into the bulk 
by a finite particle current. 

In this paper, we consider a particle non-conserving 
version of ASEP with open boundaries. Particles, that 
are injected at one boundary with a rate a, hop forward 
with rate unity (no backward hopping is allowed) with 
hard-core exclusion till they reach the other end of the 
lattice where particles are withdrawn at a rate 1 — 7 
In addition, there is particle non-conservation at the bulk 
due to the possibility of attachment (detachment) of the 
particles to (from) the lattice with rate bj a (LOd)- With 
the variation of the injection and withdrawal rates at the 
boundaries, such systems exhibit different phases which 
are characterized by the shape of the density profiles and 
the nature of the current densities^- A mean-field anal- 
ysis, exact for the hard-core case, leads to a rich phase 
diagram (see figm-c^and^ compared to systems without 
particle adsorption/ desorption kinetics (Langmuir Kinet- 
ics). For example, there is a phase where the density 
profile exhibits a jump from a low to a high value at 
some point in the bulk. This new phase known as the 
shock phase(S) appears in addition to the high(H) and 
low-density(L) phases where the bulk-density, though not 
constant, remains above and below half respectively. Al- 



though the shock phase is present for both K = 1 and 
K 7^ 1, with K = uj a /u!d, these two cases are distinctly 
different since in the first case a special particle-hole sym- 
metry allows a maximal current (M) phase with bulk den- 
sity equal to 1/2. M phase and the coexistence of M phase 
with L and/or H phase (LMH, LM and HM phases) for 
K = 1, differences in the shapes of the phase boundaries 
and in the nature of the phase transitions in the two cases 
make the problem altogether nontrivial. 

In a recent work 0] > Mukherji and Bhattachrjee 
have studied the phase-transition between low-density 
and shock phases using boundary layer techniques. Al- 
though the phase diagram has been studied numerically 
or through a direct solution of the mean-field equation 
earlier, the boundary layer analysis in [Toj | reveals a new 
feature that the transition to the S phase from the L 
phase has a precursor of a critical deconfincmcnt of a 
boundary layer near the open end. In addition, this 
approach provides a general framework for character- 
izing the transitions for different values of K and also 
for studying different mean-field models with additional 
inter-particle interactions. The purpose of the present 
work is to employ the boundary layer theory to the mean 
field version of the model to understand various phase 
transitions from a boundary layer point of view. From 
this analysis, it becomes clear that the divergences of var- 
ious length scales associated with the boundary layers in- 
dicate the emergence of bulk phase transitions. Our anal- 
ysis also reveals possibilities of surface transitions within 
the low- and high-density phases for both K = 1 and 
K 1. Except for the analysis of the surface transition, 
the discussion is primarily limited to K — 1 since, for 
K = 1, various phases and phase transitions give more 
scope for studying the effect of bulk phase transitions on 
the boundary layers. 

The paper is organized as follows. The model and its 
symmetry properties have been discussed in section II. 
We discuss the boundary layer analysis for this model in 
section III. This section is divided into three subsections 
for the discussion of low-density, high-density, and maxi- 
mal current related phases. Observations from the above 
analysis related to the behavior of the boundary layers 
near the phase boundaries are mentioned in section IV. 
Finally, we conclude with a summary of our results in 
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section V. 




FIG. 1: Phase diagram for K = ui a /uJd — 1 and Q = uidN < 
.5. iV is the number of lattice points. The coexistence of a 
maximal-current and a low- or a high-density phase is repre- 
sented by LM or HM respectively. The coexistence of low- 
density, maximal-current and high-density phases is similarly 
represented by LMH. Approach to different phase boundaries 
are indicated by the paths with arrows. "1" and "r" indicate 
presence of boundary layers at left and right boundaries re- 
spectively. Dashed lines represent surface transition lines. 
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FIG. 2: Phase diagram for K 1 and 0. = u d N < .5. N 
represents the number of lattice points. Dashed lines repre- 
sent surface transition lines, "bl" implies boundary layer. "1" 
and "r" indicate presence of boundary layers at left and right 
boundaries respectively. The filled black circle represents the 
critical point. 



II. MODEL 

To describe the ASEP of noninteracting single species 
of particles, we consider an one dimensional chain of N 
lattice points and length Denoting the occupancy of 
the ith site by r t , which assumes values or 1 depending 
on whether the site is empty or occupied by a particle, 



one may write down the mean-field master equation de- 
scribing the time evolution of ?ij = (tj), with < ... > 
denoting the statistical average, as 



drii 
~dt 



- Hi) - 71,(1 - n i+1 ) + 

UJ a (l - n.j) - LUdTlt. 



(1) 



In this equation, the mean-field approximation is imple- 
mented by neglecting correlation as 



(nn+i) = (n)(T i+ i). 



(2) 



In the large N limit, with the lattice spacing l/N — > 
0, one can go over to the continuum by substituting < 

n±i >= p(x,t) ± i|£ + f^r-... where p(x,t), is an 
average density at position x = i l/N. The choice / = 1 
simplifies the notation and restricts the variable x within 
a range {0, 1}. Keeping terms up to 0(N~ 2 ), one obtains 
the following equation describing the shape of the density 
profile in the steady state 



dx dx 



n/ (p) = o, 



(3) 



with ft = u d N, e = 1/(2JV) and f l (p) for i = 0, 1, 2 given 
as follows. For the dynamics that we have considered 
here, 

h(p) = l, h(p) = 2p-l, f (p) = K(l-p)-p. (4) 

These / functions, in general, contain information about 
the dynamics. Such general form for equation © will 
be useful laetr to make certain general predictions about 
the boundary layers. 

The last term in equation J3J originates from particle 
adsorption/desorption kinetics and is responsible for the 
loss of particle number conservation in the bulk. In the 
absence of this term, the full equation describing the time 
evolution of p(x, t), is expressible as a continuity equation 



dp(x, t) 
dt 



dj_ 

dx ' 



(5) 



with the particle current-density j given as 

dfi(p) 



dp 

j = - e fo(p)-fa - h(p)> where " J Qp J = MP)- ( 6 ) 



In case of equation the current is given by 
J = -e d £ + p{l-p). 



(7) 



In the coninuum limit (e — ► 0), the current-density j = 
p(l — p) is bounded, j < 1/4. Since, in the maximal 
current phase, the bulk density is p = 1/2, the current- 
density acquires its maximum value j — 1/4 in this phase. 
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A. symmetries 

A better understanding of the phase diagram can be 
gained from the particle-hole symmetry of the problem. 
The hopping of a particle in the forward direction is 
equivalent to the hopping of a hole in the backward direc- 
tion. Similarly, the injection of particles at one bound- 
ary with a certain rate is equivalent to the withdrawal of 
holes with the same rate. The attachment or detachment 
of particles can be interpreted as detachment or attach- 
ment of holes. The invariance of equation J2J along with 
f-functions in equation 0} under the transformation 

p — > 1 — p, x — > 1 — x (8) 
uj a ^ tod a <-> 1 — 7. (9) 

implies that the particle-hole symmetry is respected by 
the system. This symmetry is not necessarily an obvious 
property of the system and can be easily destroyed by 
additional symmetry breaking interaction terms 0, Hol | . 

The situation with K = 1 is somewhat special. If the 
adsorption/desorption kinetics were the only dynamics, 
the system would have settled in a steady-state density, 
known as the Langmuir density 

Pl =K/{K+l) (10) 

determined from fo(p) = 0. If the hopping rules respect 
conservation, then in the steady state there should be a 
homogeneous current, if this is the only dynamics. The 
corresponding density p c is determined from the zero of 
/i(p), i-e. 

Jl(Pc) = 0. (11) 

In case of particle-hole symmetry, we expect p c = 1 — p c , 
i.e., p c = 1/2 to be special. This as noted earlier is 
the maximal current state. If p c = pi, then p = p c be- 
comes a particular solution of the equation. p c is the 
density at which bulk may allow nonanalytic behavior in 
the density. This feature is useful for shock formation 
|l0j though the discontinuity is rounded by e dependent 
term in equation J2J. The adsorption/desorption dynam- 
ics need not respect this symmetry of hopping and hence 
the two densities need not be equal. K — 1 is a spe- 
cial case where the two densities become equal and the 
bulk dynamics is symmetric under the transformation 
p(x) -> 1 -p(l - x). 

III. BOUNDARY LAYER ANALYSIS 

By changing the boundary values a and 7, one may 
map out all possible steady state configurations, thereby 
constituting a "phase diagram" of the non-equilibrium 
system in the N — > 00 limit. To understand all these 
phases, we use a leading order boundary layer analysis 
that provides a systematic way to generate an uniform 
approximation of the solution of equation P| . In all the 



phases, the density profile over almost the entire space, is 
described by the solution of the first order equation ob- 
tained by ignoring the second derivative term (e — ► 0) in 
equation This solution, known as the outer solution, 
is not, in general, expected to satisfy both the boundary 
conditions. In order to satisfy the boundary conditions 
appropriately, there appear special regions with bound- 
ary layers or shocks. Description of these special regions 
requires going beyond the first order equation. The solu- 
tions describing the boundary layers or shocks are known 
as inner solutions. The constants in different solutions of 
the differential equations are determined either by the 
boundary conditions or by smooth joining of the inner 
and outer solutions. 

The low and high-density phases are related to each 
other due to particle-hole symmetry and does not require 
any separate treatment. The difference between the low- 
density and the maximal current phase arises from the 
outer solution itself. The outer solutions, solutions of the 
first order equation (equation (J3J) in the limit e — > 0), are 

pi.out = 1/2, P2,out = + c, for K = 1 (12) 
and 0.x = g(p ut) — c, with (13) 

g(j>) = yTk {2p + §TT log[/ ' ~ (1 + K)p]) (14) 

for K±\. 

Here, c is an unknown constant to be determined from 
the boundary condition. In the low- and high-density 
phases, the density profile, over almost the entire space, 
is described by the linear solution in equation l|12|) for 
K = 1 or by the solution in equation (f 1 4|) for K ^ 1. 
The maximal current phase, present only for K = 1, has 
a constant density profile described by the outer solution 
Pi, out- In the phase diagram for K = 1, there are other 
regions with coexistence of the maximal current phase 
with low- or high-density phases or both. The density 
profile in these phases has constant part pi j0 ut(%) = 1/2 
as well as linear parts p2,out = Qx + c with different parts 
joined smoothly through specific inner solutions. 

The scheme to find the inner solution varies depending 
upon the kind of matching conditions the inner solution 
has to satisfy. The procedure to find out the inner so- 
lution for L/H phases is different from that for phases 
involving M phase. The difference arises because in the 
later case, the inner solutions are required to saturate to 
p(x) = 1/2 at either of the two sides and at this value of p, 
fi(p) = 0. In the low- or high-density phases, the bound- 
ary layer is not required to saturate to 1/2. Because of 
these differences, we discuss , L, H and M-related phases 
in different subsections. 



A. Low-density phase: Surface transition 

In the low-density phase, c = a for K = 1 and c = 
17(a) for K 7^ 1 since the outer solution satisfies the left 
boundary condition. To find the inner solution, one needs 
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to express equation Q in terms of x = (x — Xd)/e, where 
Xd denotes the location of the solution. In the e — > 
limit, the inner solution is the solution of 



dx 



F( Pin ) 



where 



F{p) = f 1 { Po )-f l {p). 



(15) 



(16) 



The matching condition p- m {x — > — oo) = p = p2,out(x = 
1) for smooth joining of two solutions has already been 
incorporated. The inner region obeys the particle conser- 
vation condition and the continuity equation demands a 
homogeneous current. Consequently, the current has to 
be equal to the bulk current entering the region. This is 
the content of equations (|T5)l and ltl()|) though derived in 
a different way. 

The condition for saturation of the inner solution as 



co requires 



F(p) = for p = p s > p . 



(17) 



If 7 > p ln (x — > oo), the boundary condition cannot be 
satisfied by p m . As a result, the surface layer deconfines 
from the surface and enters into the bulk with the outer 
solution again appearing at the right edge to satisfy the 
right boundary condition. This mechanism leads to the 
formation of a shock in the density profile at the bulk. 
The bulk transition to the shock phase, therefore, has a 
precursor of deconfinement of the surface layer. Since p 
is a function of a, one has a phase boundary on the a — 7 
plane 7 = p s (p (a)) that separates the low-density phase 
from the shock phase. From equations l@J and (|16[) . we 
find p s = 1 — p a . Assuming simple zeros for F(p), we 
write 



F{p) = ~{P - Po){p - Ps)<l>{p), 



(18) 



where </>(p) can be determined using equation (jSJ). This 
form of F (p) is convenient as it yields the length scale 
associated with the crossover of the surface profile to the 
bulk profile quite generally. The large x behavior of the 
inner solution can be found from 



dx 



(P~ Ps) 
w(a) 



where 



w{a) = (p s - p ) 



-l /2 (Ps) 
<P(Ps) ' 



(20) 



Equation (|19|) shows w(a) as the characteristic length 
scale for the approach to saturation or to the bulk density. 
w(a) can be made to diverge by changing a and this 
locates a critical point on the phase boundary between 
the low-density and the shock phase. At this critical 
point (a c ,7 c ), p s = p Q = j c . 

Given this form of F(p), in the low-density phase, the 
slope of the inner solution at x = 1, i.e. at p = 7, is 



positive if p < 7 < p s and negative if 7 < p . Thus, 
with 7 SC = p 01 there is an increasing solution at x = 1 for 
7 > 7 SC and a decreasing one for 7 < 7 SC . This leads to 
a general conclusion, that for every a, if there is a bulk 
transition at 7 = p s {p {a)), there is a "boundary" tran- 
sition at 7 = p (a). This defines a surface transition line 
7 = 7 S urf(aO in the low-density phase. The transition is 
strictly at the boundary because as one crosses the phase 
boundary, the bulk density profile remains same but the 
boundary profile changes drastically. Such surface tran- 
sition is expected to be true generally for all K and also 
for interacting system, wherever a transition to a shock 
phase takes place. 

In terms of w specified above, the inner solution, in 
general, has a form 



P\n{x) = p S in (x/2w + 0, 



(21) 



with S'in — > 1 as x — > —00. The constant, £, can be 
determined from the constraint S- m (x = 0) = 7 and the 
explicit form of F(p) specified in equation (|18f) . Near the 
surface transition line £ has a logarithmic divergence 

£~]n| 7 -A,|, (22) 
which leads us to define an exponent £ s as 

£ -I Aa | Cs , for a = a surf ± Aa. (23) 
In our case, 

Cs = O(log). (24) 

In the following subsection, we illustrate this general 
approach for special cases of K = 1 and K ^ 1. 



1 . Results for K = 1 and K ^ 1 

Following the above arguments, surface transition lines 
can be obtained for all values of K . For all K , the explicit 
inner solution with positive slope at x — 1 is 

Pin(i) = 1+ (1 ~ 2 2Po) tanh[|- + £], (25) 



f\g\ where £ is a constant and 



w = 1/(1 - 2 Po ). 



(26) 



This solution appears for 7 > p a . 

As x — ► 00, the inner solution saturates to p s = 1 — 
p a . The bulk transition to the shock phase occurs when 
the saturation valuei, p s , of the surface layer, is smaller 
than 7. The phase boundary between the low-density 
and the shock-phase is, therefore, given by the equation 
1 — p {ot) = 7- The surface transition to a boundary 
layer of negative slope occurs when 7 < 7 SC = p Q . The 
boundary layer, in this case is 



1 1 — %Po , r/ x ^. 

^ + ^7^coth — + , 

2 Z Zw 



(27) 
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where w is same as in equation (|26fl . w and £ together 
determine the position of the "virtual origin" at which 
the argument of the inner solution vanishes. This is sim- 
ply in the sense of mathematical continuation since the 
virtual origin may lie well beyond the physical range of 
x, {0, 1}, with an unphysical value of the density. 




X 



FIG. 3: Density profile with positive slope at the right bound- 
ary for K = 1 with a = .1, 7 = .44, e = .0035 and ft = .2. 
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FIG. 4: Density profile with negative slope at the right bound- 
ary for K — 1 with a = .2, 7 = .3, e = .0035 and ft — .2. 



This intersection point is, therefore, also the critical point 
(a c ,7 c ) at which w(a) diverges. For K = 1, the criti- 
cal point is (1/2 — 0, 1/2). For 7 < 7 C , this boundary 
layer with negative slope at x = 1 cannot shocken if a 
is increased. However, as a is increased, this decaying 
boundary layer continues to be there in the low-density 
maximal-current (LM) phase. 

For K 7^ 1, p is the solution of g(p ) = O + g( a ) with 
g(p) given in equation 1|14|) . In this case, therefore, the 
surface transition line 7 = 7 SC = p is determined from 
the solution for 7 from the equation 

fl(7) = fi + ff(a). (32) 

The condition for shock formation from a low-density 
phase with tanh type boundary layer, on the other hand, 
leads to a phase boundary given by the solution of the 
equation 

g(l-j) = n + g(a). (33) 
The solutions for 7 from equations (|3*2^l and (|3^|l arc 
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FIG. 5: Solid line represents the surface transition line in the 
low-density phase for K — 3 and ft — .1. The dotted line 
corresponds to the phase boundary between the low-density 
and the shock phase. The intersection of the two lines is the 
critical point (a c ,j c ). 



For K = 1, using p = Q + a, we find 

w = 1/(1 -20 -2a). (28) 
Further, the constraint p ln (x = 0) = 7, leads to 

g^log L 7 -"- ' J - (29) 

1 1 — \l — a — 7 

for tanh type surface layer (see figure |3J) and 

e = lio g[ 7 - (n + a) j (30) 

2 7 — (1 — \l — a) 

for coth boundary layer (see figure 0. The surface tran- 
sition from the tanh type surface layer to the coth type 
boundary layer takes place if 7 < 7 SC = p a = Q + a lead- 
ing to a linear surface transition line, 7 = O + a, on the 
a — 7 plane. As this line is approached from either of the 
two low-density phases, £ diverges logarithmically as 



£ ~ log I 7 — O — a 



(31) 



The surface transition line p {ct) — 7 and the shock phase 
boundary p s = 1 — p Q {ct) = 7 intersect at p Q = 1/2. 



symmetric around 7=1/2 (see figure and the inter- 
section of the two solutions at 7 = j c : = 1/2 is the critical 
point (a c , 7 C ) where 



g(lc) = Q + g{a c ). 



(34) 



Since the surface transition line is symmetric to the phase 
boundary between the shock and the low-density phases, 
shapes of these to lines close the critical point are same 
jToj . The shape of the surface transition line can be ob- 
tained independently by substituting 7 = 7 C — A7 and 
a = a c — Aa in equation and expanding it in small 
A7 and Aa. This leads to a general equation 

ff '(7 c )A 7 - 5 "( 7c )(A 7 ) 2 /2 + .... = g'(a c )Aa. (35) 

Since for K ^ 1, g'{^ c ) = 0, the shape of the surface 
transition line near the critical point is given by A7 ~ 
Aa*- with 



X- 



1/2. 



(36) 



Similar to the K = 1 case, the coth boundary layer, in 
the low-density phase, leads way to a decaying boundary 
layer in the shock phase for 7 < j c . 
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The mechanism for the formation of the shock for 
7 > 7c, is different from that for 7 < 7 C , since in the 
latter case, there is no tanh type boundary layer to be 
deconfined to form a shock. The effective boundary con- 
dition for shock formation, for 7 < 7 C is same as 7 = 7 C 
since the right branch of the outer solution satisfies this 
effective boundary condition at its right edge. The orig- 
inal boundary condition p(x = 1) = 7 with 7 < j c is 
satisfied finally with a decaying boundary layer. As a 
consequence of this, a c continues to be the critical value 
of a for shock formation for all 7 < 7 C leading to a vertical 
phase boundary between the low-density and the shock 
phase. As a increases in the shock phase, the discontinu- 
ity, formed at x = 1, moves towards x = till it reaches 
the other end at the high-density-shock phase boundary. 
Whereas for 7 > 7 C , the deconfined boundary layer or 
the shock at x = 1 has a finite height, for 7 < -f c , the 
shock height increases continuously from zero to a finite 
value as a is increased. Although the mechanism of shock 
formation is different for 7 > j c and 7 < 7 C , the discon- 
tinuity is always described by a tanh type inner solution. 
Alternatively, for all values of 7, the emergence of shock, 
as one approaches the shock phase from the high-density 
side, is through the dcconfincmcnt of a tanh boundary 
layer of finite height at x = 0. As a consequence of this, 
the high-density-shock phase-boundary cannot have any 
critical point on it. 

The entire process of surface transition and then shock 
formation by changing 7 for a given a can be understood 
on a more physical ground. In the low-density phase, for 
small 7, the withdrawal rate at x = 1 is high. Since the 
bulk dynamics is completely controled by hopping and 
adsorption/desorption kinetics, a large withdrwal rate 
with a fixed a, causes particle depletion at x = 1 to 
be described by a "virtual origin" somewhere at x > 1. 
As 7 is increased, the withdrawal rate decreases and we 
reach a situation where there is neither any depletion nor 
any accummulation of particles at the end. The "virtual 
origin" is shifted to 00 now. If 7 is increased further, 
the withdrawal rate is too slow to get rid of particles 
that bulk dynamics can feed. This leads to an accum- 
mulated region at the boundary until a stage where this 
accummulatcd region becomes macroscopic. Beyond this 
point, the withdrawal rate controls the density near the 
end and the injection rate controls the remaining part of 
the density with the two parts joined through a shock or 
dicontinuity. 



B. High-density phase 

1. K = 1 

For K = 1 , the H-phase can be completely understood 
from the knowledge in L-phase by exploiting the particle- 
hole symmetry. In this sense, the high-density phase for 
particles is equivalent to the low-density phase for holes 
with density profile same as that in figure [3] or figure 0] 



when the coordinates are transformed as x — > 1 — x. In 
this transformed coordinates, above the line 7 = fi + a, 
we have the density profile for holes same as figure 
The boundary layer here is described by equation 125|) 
with the constants 

w = 1/(27 -2fi- 1) and (37) 

£ = 5log[( 7 -n-a)/(7 + a-l-n)]. ( 38 ) 
It is straightforward to verify that in original coordi- 
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FIG. 6: Density profile with positive slope at the left bound- 
ary for K = 1 with a = .55, 7 = .8, e = .001 and Q, = .2. 

nates, the particle density profile appears as figure It 
consists of a linear profile satisfying the right boundary 
condition and a boundary layer as in equation <|25[) with 

£ = ^bg[( 7 + a-l-JJ)/(7 -«-«)]. (39) 

and w same as in <|37[1 . The particle- hole symmetry also 
guarantees a surface phase transition in the high-density 
phase across 7 = Q + a line. The tanh type boundary 
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FIG. 7: Density profile with negative slope at the left bound- 
ary for K = 1 with a = .6, 7 = .76, e = .001 and Q — .2. 

layer at x = with positive slope changes to a bound- 
ary layer with negative slope (see figure UJ described by 
equation 1|27|) with constants that can be determined us- 
ing the symmetry. 

2. K^l 

The high-density phase can be divided into two major 
parts. For 7 < "f c — 1/2, the density profile has bound- 
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ary layers on both the ends. In addition to a tanh type 
boundary layer at x = 0, there is a decaying boundary 
layer at x = 1. The right boundary layer helps the den- 
sity profile satisfy the boundary condition at x = 1 from 
a value p = 1/2. Thus as in the case of shock phase 
with 7 < 1/2, the effective boundary condition on the 
right edge for the outer solution continues to be p = 1/2 
in this part of the H phase. In the other part of the H 
phase, the density profile has only one tanh boundary 
layer at x = 0. 

The transition to the shock phase from the high- 
density phase takes place through the deconfinement of 
the inner solution 

,_ 1/2+ Se^t^MziH+a ,40) 

at x = 0. Here p' Q = p2,out{x = 0) is the value of the 
outer soultion at x = 0. For 7 > 1/2, the outer solution 
obeying the boundary condition p(x = 1) = 7 can be 
obtained by solving 



C. Boundary layers in LM and HM phases for 

K = 1 

For K = 1, a = 1/2 - fi and 7 = 1/2 + Q arc the 
boundaries for the low-density and high-density phases 
respectively since tanh or coth type boundary layers are 
no more valid on these lines. The LM phase and the 
symmetrically opposite HM phase appear in the regimes 
{1/2 - n < a < 1/2,7 < I/ 2 } and {a > 1/2, 1/2 < 
7 < 1/2 + 0} respectively. In the LM phase (see figure 
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g{p' ) = -n + g( 7 ). (41) 

For 7 < 1/2, the outer solution satisfies the effective 
boundary condition p(x = 1) = 1/2 and it is the solution 
of the equation 

g(p' ) = -a + g(l/2). (42) 

The boundary layer deconfines, whenever a is smaller 
than 1 — p' a , the saturation value of the inner solution 
in equation l|40(l . This leads to the high-density-shock 
phase boundaries, 

5 (1 -a) = -Q + g(l/2) for 7 < 1/2 (43) 
g{l - a) = -Q + 5(7) for 7 > 1/2. (44) 

Since for 7 < 1/2, the value of the critical a does not 
depend on 7, the phase boundary is vertical for all 7 < 
1/2, with a K dependent value of a. 

Across the surface transition line, the slope of the 
boundary layer at x — changes sign. If the value of 
p' Q is larger than a, a boundary layer with a positive 
slope at x = is expected. In the reverse situation, one 
expects a boundary layer with a negative slope at x = 0. 
The transition lines arc, therefore, given by 

g(a) =-fl + g(l/2) for 7 < 1/2, (45) 
g(a) = -ft + 3(7) for 7 > 1/2. (46) 

The surface transition lines in equations (|45|) and H46|) 
are represented by dashed lines 1 and 2, respectively, in 
figure |5] These two lines meet at 7 = 1/2 with a value 
of a that depends on K. The surface transitions across 
both the lines are associated with the divergence of £ and 
there is no critical point on these lines. 



FIG. 8: Density profile in the low-density maximal current 
(LM) phase with a = .45, 7 = .4, e = .002 and Q, = .2. 

IHJ, the linear profile, p(x) = tlx + a, satisfies the left 
boundary condition and continues till x c \ = (1/2 — a)/£l 
where p(x c \) = 1/2. The constant profile continues till 
the other end where a boundary layer finally satisfies the 
boundary condition. In the HM phase (see figure 0), the 
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FIG. 9: Density profile in the high-density maximal current 
(HM) phase with a = .55, 7 = .58 e = .001 and = .2. 

linear profile satisfying the right boundary condition ends 
at x c h — (1 — 27 + 2f2)/2f2 with constant density profile 
for x < x c h and a boundary layer at x — 0. To obtain 
the boundary layer near x = l(or x = 0) for LM (or HM) 
phase, it is useful to express equation in terms of 

f(x*) = (2p-l)/V~e, (47) 

with x* = (x — xq) /■*/£, where xq, as before, represents 
the center of the solution. In terms of f(x*), the new 
equation is 

2 dx* 2 2dx* K ' 
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A phase plane analysis of equation (|48|l is useful to iden- 
tify the appropriate inner solution that can satisfy the 
boundary conditions. Since the details of such analy- 
sis can be found in [TlT | or in other mathematics text- 
books, we briefly mention the basic principles and obtain 
the contour-plot (figure 110(1 in an Appendix. Denoting 



df 

dx * 



p, we have 



dp 20 — p 



(49) 



To obtain the boundary layer at x = 1 for the LM 
phase, we require the specific solution which satisfies the 
boundary conditions p(x*) = 1/2, as x* — > — oo and 
p(x*) =7 as x — > 1. These conditions are fulfilled 
by the contour marked with an arrow in figure ITU1 The 
origin (/ = p = 0) is a fixed point for the differential 
equations and a simple linearization around this fixed 
point leads to p = (20) 1 / 2 /. The approach to p = 1/2 
as x* — > — oo is, therefore, exponential as 



/ ~ exp[(20) 



V2„ 



(50) 



Away from the fixed point (/ = p = 0), the solution of 
the differential equation is f 2 ~ — 2p, implying / ~ 2/x* , 
as x* — > 0. The boundary condition p(x = 1) = 7 further 
leads to x$ = 1 — 2 -\ - As 1 -> 1, the density profile 
thus has an algebraic decay as 



p{x) 



l + 2e/(2 7 -l) 



(51) 



For the HM phase, the boundary layer at the left edge 
has to be described by the specific solution that satsi- 
fies the boundary condition at x = and approaches 
exponentially as x* — > 00. It can be shown that this so- 
lution is the same as 1(5 1[) near x — 0, with xq = 
The matching between the linear and the constant pro- 
file around x c i or x c h can also be done by choosing the 
appropriate solution from the contour plot. In the case 
of LM phase, this solution should merge to the linear one 
for x* — > —00 and, approach exponentially as x* — > 00, 
with the same length scale as in equation ((50(1 . Simi- 
lar analysis is suitable also for the maximal current (M) 
phase, that appears for a > 1/2 and 7 < 1/2. In this 
phase, the density profile is fixed at 1/2 over the entire 
lattice except for decaying boundary layers same as those 
in LM and HM phases at x = 1 and x = 0, respectively. 



7 = 1/2, the linear profile satisfies both the boundary 
conditions and there is no need of any boundary layer. 
Since coth a; ~ l/i as 1 -> 0, there is an algebraically 
decaying inner solution 

^ = ^ + ,-l + 2e/(27-l) ' ^ 

near the phase boundary a = 1/2 — for 7 < j c as 
This decay is exactly the same as that we 



1 



2e 



find in the density profile in equation 1(51(1 that describes 
the boundary layer in the LM phase (on the other side of 
the phase boundary, a = l/2 — O). From this analysis, it 
becomes clear, how the coth boundary layer transforms 
to an algebraically decaying boundary layer appropriate 
for LM phase across the phase boundary. As x — ► — 00, 
the coth inner solution in the L-phase approaches p a = 
O + a exponentially as 



p(x) = + a-(l-20- 2a)e 



2[(l-2n-2a)£/2+J] 



, (53) 



with a length scale 1/(1 — 20 — 2a) that diverges as the 
phase boundary a = 1/2 — O is appraoched. It is inter- 
esting to note that as the phase boundary is approached, 
the algebraic decay of the inner solution near x = 1 is 
same in both L and LM phases. The approach of the 
two inner solutions in L and LM phases to the bulk den- 
sity, on the other hand, is exponential with two different 
length scales. Further, x diverges as xq ~ (27 — 1) _1 
as one approaches the LM and LMH phase boundary 
along path (3). Situation is symmetric as one approaches 
the HM and H phase boundary. As the phase boundary 
a > 1/2 and 7 = l/2 + O is approached from the H 
phase along path (4), £w in the coth boundary layer at 
x = appraoches l/(2a — 1). The power-law decay of 
the boundary layer near x = as 



1 

~ 2 + ~ 



2e/(l-2a) 



(54) 



is the same as that for the boundary layer in HM phase. 
As HM or LM phases are approached from M phase, the 
boundary layer at x = 1 or at x = 0, respectively, dis- 
appears and a linear profile at the respective edge starts 
appearing. This is marked by the divergences of xq as 
(27 — l)" 1 or (2a — 1) _1 along path (5) or (6) respec- 
tively. 



IV. APPROACH TO VARIOUS PHASE 
BOUNDARIES 

Our final aim is to study the boundary layers for K = 1 
as different phase boundaries are approached. As the 
phase boundary a = l/2 — O is approached from the 
low-density side (along path 2), ui£ approaches a finite 
value 1/(27 — 1) which diverges as the special point 7 = 
1/2 is approached. Since at this special point, O + a = 



V. SUMMARY 

In summary, we have shown that various bulk phase 
transitions in particle nonconserving asymmetric simple 
exclusion processes with open boundaries are associated 
with divergences of various length scales related to the 
boundary layers of the density profile. The nature of the 
divergences are obtained by implementing the techniques 
of boundary layer analysis on the mean-field equation for 
the steady-state density profile. In addition to this, we 
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also show that the low- and high-density phases exist 
in two different surface phases with distinctly different 
surface layers in the density profile. 



APPENDIX A: PHASE-PLANE ANALYSIS FOR 
EQUATION (JUJ 



eigenvalues (2ri) 1 / 2 and — (20) 1 / 2 are i (2ft) 1 1 2 J anc ^ 

(2S1) 1 / 2 ) rcs P cc tivcly. In the following, we present 

the contour-plots in f — p plane for different boundary 
conditions. 



In terms of p and /, second order equation l (4*8jl can be 
decomposed into two coupled first order equations 



dx* 



dp 

dx* 



p and 



(Al) 
(A2) 



(/ = p = 0) is a fixed point for these equations. Lin- 
earization around this fixed point leads to the following 
matrix equation 



d 
dx* 



8p 



1 

2n o 



dp 



(A3) 



which can be solved through standard diagonaliza- 
tion scheme. The eigen vectors corresponding to the 




FIG. 10: Contour plots from equations l|Al(l and ljA2fl with 
ft = .1. The arrows on different contours indicate directions 
of increasing x* . 
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